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Parvibaculum lavamentivorans DS-1^ is the type species of the novel genus Parvibaculum in 
the novel family Rliodobiaceae (formerly Ptiyllobacteriaceae) of the order Rliizobiales of Al- 
phaproteobacteria. Strain DS-1^ is a non-pigmented, aerobic, heterotrophic bacterium and 
represents the first tier member of environmentally important bacterial communities that cata- 
lyze the complete degradation of synthetic laundry surfactants. Here we describe the features 
of this organism, together with the complete genome sequence and annotation. The 
3,914,745 bp long genome with its predicted 3,654 protein coding genes is the first com- 
pleted genome sequence of the genus Parvibaculum, and the first genome sequence of a rep- 
resentative of the family Rhodobiaceae. 



Introduction 



(DSM13023 = NCIMB13966) was isolated for its 
ability to degrade linear alkylbenzenesulfonate 
[LAS), a major laundry surfactant with a world- 
wide use of 2.5 million tons per annum [1]. Strain 
DS-IT was difficult to isolate, is difficult to culti- 
vate, and represents a novel genus in the Alpha- 
proteobacteha [2,3]. Strain DS-1 catalyzes not only 
the degradation of LAS, but also of 16 other com- 
mercially important anionic and non-ionic surfac- 
tants [hence the species name lavamentivorans = 
consuming [chemicals] used for washing [3]). The 
initial degradation as catalyzed by strain DS-It 
involves the activation and shortening of the alkyl- 
chain of the surfactant molecules, and the excre- 
tion of short-chain degradation intermediates. 
These intermediates are then completely utilized 
by other bacteria in the community [4,5]. 



Parvibaculum lavamentivorans strain DS-l^ 



P. lavamentivorans DS-It is therefore an example 
of a first tier member of a two-step process that 
mineralizes environmentally important surfac- 
tants. 



Other representatives of the novel genus Parviba- 
culum have been recently isolated. Parvibaculum 
sp. strain JP-57 was isolated from seawater [6] 
and is also difficult to cultivate [3]. Parvibaculum 
indicum sp. nov. was also isolated from seawater, 
via an enrichment culture that degraded golycyc- 
lic aromatic hydrocarbons [PAH) and crude oil [7]. 
Another Parvibaculum sp. strain was isolated from 
a PAH-degrading enrichment culture, using river 
sediment as inoculum [8]. Parvibaculum species 
were also reported in a study on marine alkane- 
degrading bacteria [9]. Parvibaculum species are 
frequently detected by cultivation-independent 
methods, predominantly in habitats or settings 
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with hydrocarbon degradation. These include a 
bacterial community on marine rocks polluted 
with diesel oil [10], a bacterial community from 
diesel-contaminated soil [11], a petroleum- 
degrading bacterial community from seawater 
[12], an oil-degrading cyanobacterial community 
[13] and biofilm communities in pipes of a district 
heating system [14]. Parvibaculum species have 
also been detected in denitrifying, linear- 
nonyl^henol [NP) degrading enrichment cultures 
from NP-polluted river sediment [15] and in 
groundwater that had been contaminated by li- 
near alkyl benzenes [LABs; non-sulfonated LAS] 
[16]. Additionally, Parvibaculum species were 
detected in biofilms that degraded ^olychlori- 
nated biphenyls (PCBs) using pristine soil as in- 
oculum [17], and in a PAH-degrading bacterial 
community from deep-sea sediment of the West 
Pacific [18]. Finally, Parvibaculum species were 
detected in an autotrophic Fe[II)-oxidizing, ni- 
trate-reducing enrichment culture [19], as well as 
in Tunisian geothermal springs [20]. The wide- 
spread occurrence of Par-vibaculum species in ha- 
bitats or settings related to hydrocarbon degrada- 
tion implies an important function and role of 
these organisms in environmental biodegradation, 
despite their attribute as being difficult to culti- 
vate in a laboratory. 

Here we present a summary classification and a 
set of features for P. lavamentivorans DS-It, to- 
gether with the description of a complete genome 
sequence and annotation. The genome sequencing 
and analysis was part of the Microbial Genome 
Program of the DOE Joint Genome Institute. 

Classification and features 

p. lavamentivorans DS-It is a Gram-negative, non- 
pigmented, very small [approx. 1.0 x 0.2 |im), 
slightly curved rod-shaped bacterium that can be 
motile by means of a polar flagellum [Figure 1, 
Table 1). Strain DS-l^ grows very slowly on com- 
plex medium [e.g. on LB- or peptone-agar plates) 
and forms pinpoint colonies only after more than 
two weeks of incubation. The organism can be 
quickly overgrown by other organisms. Larger 
colonies are obtained when the complex medium 
is supplemented with a surfactant, e.g. Tween 20 
[see DSM-medium 884 [29]) or LAS [3]. When 
cultivated in liquid culture with mineral-salts me- 
dium, strain DS-l^ grows within one week with 
the single carbon sources acetate, ethanol, or 



succinate, or alkanes, alkanols and alkanoates [Cs - 
Cie); no sugars tested were utihzed [3]. 

To allow for growth in liquid culture with most of 
the 16 different surfactants at high concentrations 
[e.g. for LAS, >1 mM; see [3].), the culture fluid 
needs to be supplemented with a solid surface, e.g. 
polyester fleece or glass fibers [2,3]. The addition- 
al solid surface is believed to support biofilm for- 
mation, especially in the early growth phase when 
the surfactant concentration is high, and the or- 
ganism grows as single, suspended cells [non- 
motile) during the later growth phase. Growth 
with a non-membrane toxic substrate [e.g. ace- 
tate) is independent of a solid surface, and consti- 
tutes suspended, single cells [motile). We presume 
that the biofilm formation by strain DS-l^ is a pro- 
tective response to the exposure to membrane- 
solubihzing agents [cf [30]). 

Based on the 16S rRNA gene sequence, strain DSl^ 
was described as the novel genus Parvibaculum, 
which was originally placed in the family Phyllo- 
bacteriaceae within the order Rhizobiales of Al- 
phaproteobacteria [3,31]. The nearest well- 
described organism to strain DS-It is Afifella ma- 
rina [formerly Rhodobium marinum) [92% 16S 
rRNA gene sequence identity), a photosynthetic 
purple, non-sulfur bacterium. The genus Rhodo- 
bium was later re-classified as a member of the 
novel family Rhodobiaceae [26,32], together with 
two novel genera of other photosynthetic purple 
non-sulfur bacteria [Afifella and Roseospirillum), 
as well as with two novel genera of heterotrophic 
aerobic bacteria, represented by the red- 
pigmented Anderseniella baltica [gen. nov., sp. 
nov.) [33,34] and non-pigmented Tepidamorphus 
gemmatus [gen. nov., sp. nov.) [35,36]. A phyloge- 
netic tree [Figure 2) was constructed with the 16S 
rRNA gene sequence of P. lavamentivorans DS-l^ 
and that of [i) other isolated Parvibaculum strains, 
[ii) representatives of other genera within the 
family Rhodobiaceae, [iii) representatives of the 
genera in the family Phyllobacteriaceae, as well as, 
[iv) representatives of other families within the 
order Rhizobiales. The phylogenetic tree shows 
now the placement of Parvibaculum species within 
the family Rhodobiaceae, and that the Parvibacu- 
lum sequences clustered as a distinct evolutionary 
lineage within this family [Figure 2). This classifi- 
cation of Parvibaculum has been adopted in the 
Ribosomal Database Project [RDP) and SILVA 
rRNA Database Project, but not in the GreenGenes 
database. The family Rhodobiaceae has also not 



http://standardsingenomics.org 



299 



Parvibaculum lavamentivorans type strain (DS-1T) 



been included in the NCBI-taxonomy, IMG- 
taxonomy, and GOLD databases. 

Currently, 360 genome sequences of members of 
the order Rhizobiales of Alphaproteobacteria have 
been made available [GOLD database; August 
2011), and within the family Phyllobacteriaceae 
there are 21 genome sequences available 
[Chelativorans sp. BNCl, Hoeflea phototrophica 
DFL-43, and 18 Mesorhizobium strains). No ge- 
nome sequences currently exist for a representa- 
tive of the novel family Rhodobiaceae, except of 
the genome of P. lavamentivorans DS-It. 

Chemotaxonomy 

Examination of the respiratory lipoquinone com- 
position of strain DS-It showed that ubiquinones 
are the sole respiratory quinones present, and the 
major hpoquinone is ubiquinone 11 (Qll) [3]. The 
fatty acids of P. lavamentivorans are straight chain 
saturated and unsaturated, as well as ester- and 
amide-linked hydroxy-fatty acids, in membrane 
fractions [3]. The major polar lipids are 



phosphatidyl glycerol, diphosphatidyl glycerol, 
phosphatidyl ethanolamine, phosphatidyl choline, 
and two, unidentified aminolipids; the presence of 
the two additional aminolipids appears to be dis- 
tinctive of the organism [3]. The G+C content of 
the DNA was determined to be 64% [3], which 
corresponds well to the G+C content observed for 
the complete genome sequence [see below). 

Genome sequencing information 

Genome project history 

The genome was selected for sequencing as part 
of the U.S. Department of Energy - Microbial Ge- 
nome Program 2006. The DNA sample was sub- 
mitted in April 2006 and the initial sequencing 
phase was completed in October 2006. The ge- 
nome finishing and assembly phase were com- 
pleted in June 2007, and presented for pubhc 
access on December 2007; a modified version was 
presented in February 2011. Table 2 presents the 
project information and its association with MIGS 
version 2.0 comphance [39]. 




Figure 1. Scanning electron micrograph of P. lavamentivorans DS-1^. Cells derived from a liquid culture 
that grew in acetate/mineral salts medium. 
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Figure 2. Phylogenetic tree of 1 6S rRNA gene sequences showing the position of P. lavamentivorans DS-1^ relative 
to other type strains within the families Rhodobiaceae, Phyllobacteriaceae and other families in the order Rhizo- 
biales (see the text). Strains within the Rhodobiaceae and Phyllobacteriaceae shown in bold have genome projects 
underway or completed. The corresponding 16S rRNA gene accession numbers (or draft genome sequence iden- 
tifiers) are indicated. The sequences were aligned using the GreenGenes NAST alignment tool [37]; neighbor- 
joining tree building and visualization involved the CLUSTAL and DENDROSCOPE software [38]. Caulobacterales 
sequences were used as outgroup. Bootstrap values >30 % are indicated; bar, 0.01 substitutions per nucleotide po- 
sition. 
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Table 1. Classification and general features of Parvlbaculum lavamentivorans DS-1^. 



MIGS ID Property Term Evidence code" 







Domain Bacteria 


TAS 


[21] 






Phylum Proteobacteria 


TAS 


[22] 






Class Alphaproteobacteria 


TAS 


[23,24] 




Current classification 


Order Rhizobiales 


TAS 


[23,25] 






Family Rhodobiaceae 


TAS 


[23,26] 






Genus Parvibaculum 


TAS 


[3] 






Species Parvibaculum lavamentivorans 


TAS 


[3] 






Type strain DS-1 








Gram stain 


negative 


TAS 


[3] 




Cell shape 


small rod 


TAS 


[3] 




Motility 


motile 


TAS 


[3] 




Sporulation 


non-sporulating 


TAS 


[3] 




Temperature range 


mesophile 


TAS 


[3] 




Optimum temperature 


30 °C 


TAS 


[3] 




Carbon source 


acetate, ethanol, pyruvate, succinate, alkanes (Cg- 
Cig), various anionic and non-ionic surfactants 


TAS 


[2,3,5,7] 




Energy source 


chemoorganotroph 


TAS 


[3] 




Terminal electron receptor 


molecular oxygen 


TAS 


[3] 


MIGS-6 


Habitat 


aerobic habitat 


TAS 


[2,27] 


MIGS-6.3 


Salinity 


0 - 3% NaCI 


TAS 


[3] 


MIGS-22 


Oxygen requirement 


aerobic 


TAS 


[3] 


MIGS-15 


Biotic relationship 


free-living 


TAS 


[3] 


MIGS-14 


Pathogenicity 


none 

isolated from a surfactant-degrading laboratory trick- 


TAS 


[3] 


MIGS-4 


Geographic location 


ling filter that was inoculated with sludge of an indus- 
trial sewage treatment plant in Ludwigshafen, 
Germany 


TAS 


[2] 




■JdlllUlC; l^WIICri^LlWlI LllllC 


1 qqq 


TAS 




MIGS-4.1 


Latitude 


49.48 


TAS 


[2] 


MIGS-4.2 


Longitude 


8.44 


TAS 


[2] 


MIGS-4.3 


Depth 








MIGS-4.4 


Altitude 


96 m 


TAS 


[2] 



a) Evidence codes - IDA: Inferred from Direct Assay; TAS: Traceable Author Statement (i.e., a direct report exists in the 
literature); NAS: Non-traceable Author Statement (i.e., not directly observed for the living, isolated sample, but based 
on a generally accepted property for the species, or anecdotal evidence). These evidence codes are from of the Gene 
Ontology project [28]. 
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Table 2. Project information 



MIGS ID 


Property 


Term 


MIGS-31 


Finishing quality 


Finished 


MIGS-28 


Libraries used 


3.5 kb, 9 kb and 37 kb DNA libraries 


MIGS-29 


Sequencing platforms 


Sanger 


MIGS-31. 2 


Sequence coverage 


16x 


MIGS-30 


Assemblers 


Phred/Ph rap/Con sed 


MIGS-32 


Gene calling method 


Glimmer/Criteria 




Genbank ID 


17639 




Genbank Date of Release 


July 31, 2007 




GOLD ID 


Gc00631 


MIGS-13 


Source material identifier 


DSM 13023 = NCIMB 13966 




Project relevance 


Biodegradation, biotechnological 



Growth conditions and DNA isolation 

P. lavamentivorans DS-It was grown on LB agar 
plates [2 weeks) and pinpoint colonies were 
transferred into selective medium (1 mM 
LAS/minimal salts medium; with glass-fiber sup- 
plement, 5-ml scale [3]). This culture was sub- 
cultivated to larger scale [100-ml and 1-liter 
scale) in 30 mM acetate/minimal salts medium; 
cell pellets were stored frozen until DNA prepara- 
tion. DNA was prepared following the JGI DNA 
Isolation Bacterial CTAB Protocol [40]. 

Genome sequencing and assembly 

The genome of P. lavamentivorans DS-It was se- 
quenced at the Joint Genome Institute [JGI) using 
a combination of 3.5 kb, 9 kb and 37 kb DNA li- 
braries. All general aspects of library construc- 
tion and sequencing performed at the JGI can be 
found at the JGI website [41]. Draft assembhes 
were based on 76,870 reads. Combined, the reads 
from all three libraries provided 16-fold coverage 
of the genome. The Phred/Phrap/Consed soft- 
ware package [42] was used for sequence assem- 
bly and quahty assessment [43-45]. After the 
shotgun stage, reads were assembled with paral- 
lel phrap (High Performance Software, LLC). 
Possible mis-assemblies were corrected with 
Dupfinisher [46], PGR amplification, or transpo- 
son bombing of bridging clones [Epicentre Bio- 
technologies, Madison, WI, USA). Gaps between 
contigs were closed by editing in Consed, custom 
primer walk or PGR amplification [Roche Applied 
Science, Indianapolis, IN, USA). A total of 24 pri- 
mer walk reactions were necessary to close gaps 
and to raise the quality of the finished sequence. 



The completed genome assembly contains 
76,885 reads, achieving an average of 16-fold 
sequence coverage per base with an error rate 
less than 5 in 100,000. 

Genome annotation 

Genes were identified using a combination of 
Critica [47] and Ghmmer [48] as part of the ge- 
nome annotation pipeline at Oak Ridge National 
Laboratory [ORNL), Oak Ridge, TN, USA, followed 
by a round of manual curation. The predicted 
CDSs were translated and used to search the Na- 
tional Center for Biotechnology Information 
[NCBI) non-redundant database, UniProt, TIGR- 
Fam, Pfam, PRIAM, KEGG, COG, and InterPro da- 
tabases; miscellaneous features were predicted 
using TMHMM [49] and signalP [50]. These data 
sources were combined to assert a product de- 
scription for each predicted protein. The tRNAS- 
canSE tool [51] was used to find tRNA genes, 
whereas ribosomal RNAs were found by using 
BLASTn against the ribosomal RNA databases. 
The RNA components of the protein secretion 
complex and the RNaseP were identified by 
searching the genome for the corresponding 
Rfam profiles using INFERNAL [52]. Additional 
gene prediction analysis and manual functional 
annotation was performed within the Integrated 
Microbial Genomes [IMG) platform [41] devel- 
oped by the Joint Genome Institute, Walnut 
Creek, CA, USA [53]. 
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Genome properties 

The genome of P. lavamentivorans DS-It compris- 
es one circular chromosome of 3,914,745 bp 
(62.33% GC content) (Figure 3), for which a total 
number of 3,714 genes were predicted. Of these 
predicted genes, 3,654 are protein-coding genes, 
and 2,723 of the protein-coding genes were as- 
signed to a putative function and the remaining 
annotated as hypothetical proteins; 18 pseudo- 
genes were also identified. A total of 60 RNA genes 
and one rRNA operon are predicted; the latter is 
reflective of the slow growth of P. lavamentivorans 
DS-IT [54,55]. Furthermore, one Clustered Regu- 
larly Interspaced Short Palindromic Repeats ele- 
ment (CRISPR) including associated protein genes 
were predicted. The properties and the statistics 
of the genome are summarized in Table 3, and the 
distribution of genes into COGs functional catego- 
ries is presented in Table 4. 

Metabolic features 

The genome of P. lavamentivorans encodes com- 
plete pathways for synthesis of all proteinogenic 



amino acids and essential co-factors, and the cen- 
tral metabolism is represented by a complete 
pathway for the citrate cycle, glycoly- 
sis/gluconeogenesis, and the non-oxidative 
branch of the pentose-phosphate pathway; no 
candidate genes for the oxidative branch of the 
pentose-phosphate pathway or for the Entner- 
Doudoroff pathway are predicted. 

P. lavamentivorans DS-V does not grow on D- 
glucose, D-fructose, maltose, D-mannitol, D- 
mannose, and N-acetylglucosamine [3,7], and 
there are no valid candidate genes predicted in the 
genome for ATP-dependent sugar uptake systems 
or for D-glucose uptake via a phosphotransferase 
system. Similarly, no valid candidate genes were 
predicted for ATP-dependent amino-acid and 
di/oligo-peptide transport systems or for other 
amino-acid/peptide transporters, which reflects 
the poor growth of strain DS-It in complex me- 
dium (LB-medium). 



Table 3. Nucleotide and gene count levels of the genome of P. lavamentivorans DS-1^ 
Attribute Value % of total' 



Genome size (bp) 


3,914,745 


100 


DNA coding region (bp) 


3,535,064 


90.30 


G+C content (bp) 


2,439,986 


62.33 


Number of replicons 


1 




Extrachromosomal elements 


0 




Total genes 


3,714 


100 


RNA genes 


60 


1.62 


rRNA operons 


1 




Protein-coding genes 


3,654 


98.38 


Pseudo genes 


18 


0.48 


Genes with function prediction 


2,723 


73.32 


Genes in paralog clusters 


620 


16.69 


Genes assigned to COGs 


2,904 


78.19 


Genes assigned to Pfam domains 


3,054 


82.23 


Genes with signal peptides 


717 


19.31 


Genes connected to KEGG pathways 


1,085 


29.21 


Genes with transporter classification 


430 


11.58 


Genes with transmembrane helices 


782 


21.06 


CRISPR count 


1 


% of total ^ 



") The total is based on either the size of the genome in base pairs or the total number 
of protein coding genes in the annotated genome. 



304 



Standards in Genomic Sciences 



Schleheck et al. 



Table 4. Number of genes associated with the general COG functional categories in P. lavamentivorans DS-1^ 



Code 


Value 


%age 


Description 


J 


163 


5.07 


Translation, ribosomal structure and biogenesis 


A 


1 


0.0 


RNA processing and modification 


K 


243 


7.0 


Transcription 


L 


137 


3.9 


Replication, recombination and repair 


B 


1 


0.0 


Chromatin structure and dynamics 


D 


25 


0.7 


Cell cycle control, mitosis and meiosis 


Y 


0 


0.0 


Nuclear structure 


Z 


0 


0 


Cytoskeleton 


w 


0 


0 


Extracellular structures 


V 


85 


2.4 


Defense mechanisms 


T 


118 


3.4 


Signal transduction mechanisms 


M 


131 


3.8 


Cell wall/membrane biogenesis 


N 


6 


0.2 


Cell motility 


U 


39 


1.1 


Intracellular trafficking and secretion 


o 


77 


2.2 


Posttranslational modification, protein turnover, chaperones 


c 


153 


4.4 


Energy production and conversion 


G 


294 


8.4 


Carbohydrate transport and metabolism 


E 


214 


6.1 


Amino acid transport and metabolism 


F 


79 


2.3 


Nucleotide transport and metabolism 


H 


110 


3.2 


Coenzyme transport and metabolism 


1 


73 


2.1 


Lipid transport and metabolism 


P 


152 


4.4 


Inorganic ion transport and metabolism 


Q 


30 


0.9 


Secondary metabolites biosynthesis, transport and catabolism 


R 


318 


9.1 


General function prediction only 


S 


200 


5.7 


Function unknown 




1082 


31.0 


Not in COGs 



For the assimilation of acetyl-CoA from the degra- 
dation of alkanes and surfactants [2,3,5], or during 
growth with acetate, the genome of P. lavamenti- 
vorans encodes the glyoxylate cycle [isocitrate 
lyase, Plav_0592; malate synthase, Plav_0593) to 
generate succinate for the synthesis of carbohy- 
drates. The genome also encodes the complete 
ethyl-malonyl-CoA pathway to assimilate acetate 
[56]. This observation, i.e. glyoxylate cycle and 
ethyl-malonyl-CoA pathway in the same organism, 
has been made before [57], and these two path- 
ways in P. lavamentivorans DS-l^ might be diffe- 
rentially expressed under varying environmental 
conditions. 

For the degradation of alkanes and surfactants 
through abstraction of acetyl-CoA [54], the ge- 
nome contains a wealth of candidate genes for the 
entry into alkyl-chain degradation [omega- 



oxygenation to activate the chain) supplemented 
by a variety of genes predicted for omega- 
oxidations [to generate the corresponding fatty- 
acids) and fatty-acid beta-oxidations [to excise 
acetyl-CoA units). We are currently exploring this 
high abundance of genes for alkane/alkyl- 
utilization in strain DS-It by transcriptional and 
translational analysis [unpublished]. For example, 
at least nine cytochrome-P450 [CYP) alkane mo- 
nooxygenase [C0G2124), 44 alcohol dehydroge- 
nase [COG1028), 11 aldehyde dehydrogenase 
[COG1012), 20 acyl-CoA synthetase [COG0318), 
40 acyl-CoA dehydrogenase [COG1960), 31 enoyl- 
CoA hydratase [COG1024), 14 acyl-CoA acetyl- 
transferase [COG0183), six thioesterase 
[COG0824), and 17 putative long-chain acyl-CoA 
thioester hydrolase [PF03061) candidate genes 
are predicted in the genome. 
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Figure 3. Graphical circular map of the genome of P. lavamentivorans DS-1^. From outside to center; Genes on 
forward strand (color by COG categories), genes on reverse strand (color by COG categories), RNA genes (tRNA, 
green; rRNA, red; other RNAs, black), GC content, GC skew. 



Other predicted oxygenase genes comprise three 
putative Baeyer-Vilhger-type FAD-binding mo- 
nooxygenase genes [COG2072). Cyclohexanone 
and hydroxyacetophenone, which are putative 
substrates for such oxygenases (e.g [58,59]) were 
tested as carbon source for growth of strain DS-It 
as well as cycloalkanes [Ce, Cs, C12), however, none 
supported growth. The terpenoids camphor [for 
the involvement of a cytochrome- P450 oxygenase 



in the degradation pathway [60]) and geraniol, 
citronellol, linalool, menthol and eucalyptol [for 
the involvement of acyl-CoA interconversion en- 
zymes in the degradation pathways) as substrates 
for growth were also tested negative. 

In contrast to the high abundance of genes for ali- 
phatic-hydrocarbon degradation, the genome con- 
tains few genes for aromatic-hydrocarbon degrada- 
tion. One gene set for an aromatic-ring dioxygenase 
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component [Plav_1761 and 1762; BenAB-type), 
three aromatic-ring monooxygenase component 
genes (Plav_1541 and 0131, MhpA-type; Plav_1785, 
HpaB-type), and three vahd candidate genes for 
extradiol ring-cleavage dioxygenase (Plav_1539 [61] 
and 1787, BphC-type; Plav_0983, LigB-type) were 
predicted in the genome. Strain DS-l^ did not grow 
with benzoate, protocatechuate, phenylacetate. 
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